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We consider a 5D BPS dilatonic two brane model which reduces to the Randall-Sundrum model or 
the Hoi'ava-Witten theory for a particular choice of parameters. Recently new dynamical solutions 
were found by Chen et at, which describe a moduli instability of the warped geometry. Using a 
4D effective theory derived by solving the 5D equations of motion, based on the gradient expansion 
method, we show that the exact solution of Chen et al. can be reproduced within the 4D effective 
theory and we identify the origin of the moduli instability. We revisit the gradient expansion method 
with a new metric ansatz to clarify why the 4D effective theory solution can be lifted back to an 
exact 5D solution. Finally we argue against a recent claim that the 4D effective theory allows a 
much wider class of solutions than the 5D theory and provide a way to lift solutions in the 4D 
\^ • effective theory to 5D solutions perturbatively in terms of small velocities of the branes. 

o ■ 
o 

CN ■ I. INTRODUCTION 

^ . . 

^ ' The idea of using the degrees of freedom of extra spatial dimensions to unify gravity and the other interactions 
has attracted much interest. Especially, M-Theory has the potential to achieve this ultimate goal. Hofava and 
Witten showed that M-Theory compactified as /Z2 x M^" reduces to a known string theory {Es x Es Heterotic 
String) Q . From a low energy effective theory point of view, this model has a 5D spacetime (the bulk) with two 4D 
boundary hypersurfaces (the branes) Q. The remaining six spatial dimensions are assumed to be compactified. All 
the Standard Model particles are confined to the branes while gravity can propagate in the bulk. As a consequence 
QQ • of the compactification of the six spatial dimensions, a 5D effective scalar field appears in the bulk, which describes 
\^ ' the volume of the compactified 6D space. 

A novel aspect of this model is that the fifth dimension is not homogeneous. The line element for a static solution 
CN ' is given in the form dsl = h°'{r)dr'^ + h^{r)dsl{x), where r is the fifth coordinate, a, (3 are constants and h is called 
i the warp factor. This warped geometry is used by Randall and Sundrum Q to address the hierarchy problem. They 
have also showed that it is possible to localize gravity around the brane with the warped geometry, providing in this 
way an alternative to compactification 4] (see Ref. |5| for a review). 
• Recently, a new class of dynamical solutions that describes an instability of the warped geometry has been found. 
I Chen et al. noticed that it is possible to obtain a dynamical solution by replacing the constant modulus in the warp 
O-i, factor /i by a linear function of the 4D coordinates ^ . This solution describes an instability of the model as the brane 
^ ' will hit the singularity in the bulk, where /i = 0. This kind of solution exists also in the lOD type IIB supergravity 

^ ' In this paper, we study the moduli instability in a dilatonic two brane model, where the potentials for the scalar 
k>( field on the brane and in the bulk obey the BPS condition i^j. For particular values of parameters we retrieve either 
■ the Hofava- Witten theory or the Randall-Sundrum model. We identify the origin of the moduli instability using a 4D 
d I effective theory derived in Ref. (see Refs. and for different approaches). This effective theory is derived 
by solving the 5D equations of motion using the gradient expansion method , where we assume that the velocities 
of the branes are small compared with the curvature scale of the bulk determined by the warp factor. 

Despite the fact that the 4D effective theory is based on a slow-motion approximation, we will show that the 5D 
exact solutions can be reproduced in the 4D effective theory. In order to understand the relation between 4D solutions 
in an effective theory and full 5D solutions, we revisit the gradient expansion method by employing a new metric 
ansatz. Using this metric ansatz, we can clearly see why the moduli instability solution in the 4D effective theory can 
be lifted to an exact 5D solution. 

We also comment on a recent claim that the 4D effective theory allows a much wider class of solutions than the 5D 
theory fl3\. We disagree with that conclusion and we show that it is based on the restricted form of the 5D metric 
ansatz used in Ref. [l^- Using our metric ansatz, we provide a way to lift solutions in the 4D effective theory to 5D 
solutions perturbatively in terms of small velocities of the branes. 

The structure of the paper is as follows. In section II, the model used in this paper is described in detail. In section 
III, we identify the solution in the 4D effective theory that describes the moduli instability. In section IV, we revisit 
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the gradient expansion method to derive the 4D effective theory. We propose a new metric ansatz which is useful to 
relate 4D solutions in the effective theory to 5D solutions. Using this formalism we explain why the 4D solution for 
the moduli instability can be lifted to an exact 5D solution. In section V, we comment on the arguments against the 
4D effective theory. Section VI is devoted to conclusions. 

II. THE MODEL 

Our model consists of a 5D spacetime (bulk) filled with a scalar field. The fifth dimension is a compact space 
with a Z2 symmetry (i.e. identification r — r) and will be parameterized by the coordinate r, —L < r < L. The 
bulk action takes the form 



Sbulk n 



^ fd^X^^ pi? - VmV'V^V + Vbulki^)] , (1) 



where kg denotes the 5D gravitational constant. Throughout this work latin indices can have values in {r,t, x,y, z}, 
while greek indices do not include the extra dimension coordinate r. The scalar field potential will be 

^2 _ 2\ ^-2V2bv^2 



VbM^) = -(b'- e-'^^'^a', (2) 

where a and b are the remaining parameters of this model. We can retrieve either the Randall-Sundrum model (6 = 0) 
or the Hofava-Witten model (6 = 1) according to the value of b. The orbifold /Z2 has two fixed points, at r = 
and r — L, and we can put branes there. The branes' action is 

Sbranes — 5*+ + 5—, (3) 

where S+ and S*- are the positive and negative tension brane action, respectively. They are given by 

S±^T^ f d^xV^,e-^'^a. (4) 

The action of our model is 

Stotal — Sbulk "I" ^branes- (5) 

III. EXACT SOLUTIONS FOR MODULI INSTABILITY IN 4D EFFECTIVE THEORY 

A. The exact 5D solution 

Exact static solutions for the equations of motion that result from jSJ were obtained in ^8] . Recently, an interesting 
dynamical solution was found by Chen et al. 6] . They noticed that if we add a linear function of time to the warp 
factor of the solution of 8], it would still be a solution of the equations of motion. Their exact 5D solution reads 

6b^ -2 2 
„ / \r\ \ 3b^ + l „ / \r\ \ 36^ + 1 

dSi = (hr-Y) "^(''^"V) Vf^^^dx^dx", (6) 

3V2b ( \r\\ 3 V2 

where h is an arbitrary constant. 

From the above equations we can read the scale factor of the positive tension brane (r = 0) as 

ai(r) = (M^, (8) 

and the scalar field as 

'P(r- = 0,r) = |^ln(M. (9) 
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Let us choose h < and r < 0. In this case the proper distance between branes, the so-called radion, is decreasing 
according to 



dr 



1 



6&2 



(hr) 



L/l \ 3''^ + l 

757 



(10) 



When T ~oo, TZ ^ L (hr)^^ and for r = L/hl, TZ = /(Sb^ + 1) {L/l)^^ /Qb'^. Before the two branes collide, 
a curvature singularity will appear in the negative tension brane (r = i) at r = L/hl. This singularity will move 
towards the positive tension brane and will reach it at r = 0. This event represents the total annihilation of the 
spacetime. 

It is useful to note that if we drop the modulus sign in the bulk spacetime is a static black brane or black hole 
solution depending on h, see for more details. If 6 = 1, there is a timclike curvature singularity at hr = \r\/l. From 
the static bulk point of view, the two branes are moving in this static bulk and the negative tension brane first hits 
the singularity. Even if the bulk spacetime is static, the existence of the branes, which gives the modulus sign for r, 
makes the spacetime truly time dependent. 

Moduli instability is a serious problem for these types of model. In this work, we will try to understand it from a 
4D effective theory viewpoint. 



B. The exact solutions in 4D effective theory 



In 13], Kobayashi and Koyama applied the gradient expansion method to solve perturbatively the 5D equations of 
motion resulting from (|SJ|. Their 0*'* order solution reads 



dSi 



ip{r,t) = -^InJ'ir^t) + <j){t), J'{r,t) = 1 - V2b'^d{t)a\r\. 



(11) 
(12) 



The main result of the paper is a set of equations involving only four dimensional quantities that describe the dynamics 
of the unknown functions d{t), h^^{t) and (j){t), appearing in Eqs. (|ll|l and H12f) . It can be shown that their dynamical 
equations can be consistently deduced from the following action 



■ipR{h) - 



1 



2(1 + 3^2) l-^p 



(13) 



where here Vq means the covariant derivative with respect to h^^ and ^ = 1 — (1 — d) 3''^ . This is the action we 
obtain if we substitute Eqs. (|ll|l and H12|) into and integrate over the extra dimension [To| . 
Performing the conformal transformation /i^^ = e~^''"^/^i,/-0 and defining new scalar fields as 



ip = 1 — tanh 



1 + 362 
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(14) 



VI + 362 

we arrive at the action in the Einstein frame 

/ 



+ X\nip, A 



V2 362 + 1' 



2k^ 



d^x^ i?(/) _ e|„ei" - vi/|„vi;i" 



(15) 



(16) 



It is clear that the moduli fields have no potential and this is the origin of the instability. We are interested in 
cosmological solutions so we choose to be a flat FRW metric. The Einstein equations resulting from the action 
(Unj are 



« 1 /t, 

— = — U*' 
a 6 V 



2 . e'2 



-1 =i(*- + e' 



(17) 
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and the equations of motion for the scalar fields are 



*" + 2-*' = o, e" + 2-e' = o, 

a a 



(18) 



where a(r) is the FRW scale factor and the prime denotes derivative with respect to the conformal time r. Equations 
(|17I18|I can be easily integrated to find 



V oo^ + 1 



(19) 



(1 + 36^ 
6 



0^ = 1 + 36^ 



(20) 



where a, (3, 7, 5, Q and ^ are integration constants. The solutions for the fields in the original frame can be obtained 
as 



d(r) = 1 - 



{flT + 1^)2" + 1 



36^ + 1 



(21) 



In 



9, 



(22) 



where X(b) is defined in Eq. p5|) and integration constants fi, v, g are redefined from ^, 7 and 6. The radion is 
calculated as 



7^(r) 



1 



d{T)e 



\/2fc0(r) 



and the square of the scale factor on the positive tension brane is 



(23) 



(24) 



We have found the remarkable fact that for a particular choice of the integration constants we can reproduce the 
exact solution found by Chen et al. described in the previous section. If we choose integration constants obeying the 
relations 



1 R 



(25) 



we get the same 4D quantities (scale factor and scalar field on the positive tension brane) and radion as the Chen et 
al. solution, Eqs. pfel lTUIl . 
For this solution, the scale factor in the Einstein frame is given by 



a2(r) 



I \ 21 



(26) 



If we take h < 0, t < 0, this corresponds to a collapsing universe, due to the kinetic energy of the scalar fields. At 
T = L/2hl < 0, the universe in the Einstein frame reaches the Big-bang singularity. However, we should be careful to 
interpret this singularity. In fact, in the original 5D theory, r — L/2hl docs not correspond to any kind of singularity. 
In 5D theory, the negative tension brane hits the singularity at r = L/hl and the positive tension brane hits the 
singularity at r = 0. In fact at r = L/2hl, = and so the conformal transformation becomes singular and the 
Einstein frame metric loses physical meaning. 



IV. GRADIENT EXPANSION METHOD WITH A NEW METRIC ANSATZ 

In the previous section, we found that the exact solution derived in can be reproduced from the 0*'* order of 
the perturbative method. In this section, we will revisit the gradient expansion method which is used to derive the 
effective theory to clarify the reason why the exact solution derived in ^ can be reproduced within the 4D effective 
theory. 
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A. 5D Equations 



In this subsection we study the 5D equations of motion. From the action ((Sj we derive the Einstein's equations 



/.95 



g^uS^A^Be-'^'^'Sir) + ^g^J'^S^^e-^'^Sir - L) 

95 



(27) 



Assuming that the 5D hne element has the form 

dSi = gAB{X)dX^dX^ = grr{x,r)dr^ + g^j,^{x,r)dx^' dx" , (28) 

{x dependence means dependence of all the other coordinates {x, y, z, t}, except the extra dimension r) we can extract 
the junction conditions for the metric tensor from Einstein's equation, 

where the extrinsic curvature K^^, is defined by if^^ — —\drg^u- The scalar field equation of motion is 



r=0+ 



r=0+ ' 

r=L^ 



(29) 



□(^ + \/2fe ( - - ) a^e-^Vs'"^ ^ -2foa \^/^e-^^^ 5{r) ~ ^e-^^^5{r - L) 



and so the junction conditions for the scalar field are 



r=0+ 
r=L- 



^bcre 



r=0+ ■ 



In order to proceed we shall assume that the 5D line element has further symmetries, described by 



dSi 



-^^'"^'-'''''^ H^!^ {x,r)dr^ + {x,r)g^^dx^'dx'' , H{x,r) = C{x) 



I' 



(30) 



(31) 



(32) 



In this section, we will assume that the position of the second brane is r = L = L For the scalar field we will assume 
the form 



= OU2 , -I Ini? + r). 



(33) 



362 + 1 

This metric ansatz is inspired by the time dependent solution (|6I7|) of Chen et al. @. Their solution was found by 
replacing the modulus parameter in the static solution by a linear function of time. In the same manner, we introduce 
an X dependence in H through the modulus parameter C{x) in a covariant way . We also introduce the function (p 
for the scalar field moduli. In order to satisfy the junction conditions l|31|) we must have the exponential factor in the 
grr metric component. The tensor ^^i, is left completely general. 

After some mathematical manipulations of the Einstein equations we obtain 

1 



V 362 + 1 



C, 



V2h^\'l + (26^ + 1) ^l^^i. 



V26 ^-3f-3 



362 + 1 



\u ' 3^,2 4. 1 W u 

(36^ + 1) (^i^q, + + ^i"q„5^ 



(34) 



1362 



; 362 + 1 



3-27261^^ 36^ + 1 k + e^2V2b0^-j; 



k^'Pk, 



a/3 



3^ / 362 - 1 



|Q 



1 2V26 

7 362 + 1 
6^/263 



iJ^i - \/2hk 



362 _^ ]_ |a 3^2 



(35) 
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' 3V26 
362" 



(36) 



where K^j^i, is defined like K^i, — — ^9r5pi/, K is the trace of Kf_iu, \a denotes covariant derivative with respect to 5^1/ 
and R^^{g) is the Ricci tensor of g^^. Equation (|30|) transforms into 



1 1 



- \/2& ((^ ^)^ - ^ ^ii' + 7—^ -i?"^ (3V2bk + (%2 - 5) r 

i 00^ + 1 \ 



With this particular ansatz the junction conditions (|29I31|I are significantly simplified and read 



r=0+ 

r=r 



= 0, 



(37) 



(38) 



r=r 



0. 



(39) 



It is impossible to solve these equations in general, so in the next section, we will solve them using the gradient 
expansion method up to first order in the perturbations. 



B. The gradient expansion 

1. The approximation 



In order to use perturbation theory to solve a system of differential equations we need to identify the characteristic 
scale of the different terms involved in the equations and then see if there is a small parameter. 

The derivatives along the extra dimension of the conformal metric g^i, as well as the derivative of (p are of order 
l/l. We assume that variations along the branes' coordinates are small in comparison with l/l. This implies that the 
radion changes slowly or that the speed of the branes is small. More precisely, our small parameters will be 



l^dX . .) < 1 and (W, ...)'< 1, 

where . . . represents the conformal metric functions, C or (p. 

As in the usual perturbation method, we expand the unknown functions in a series 

(0) (1) 

g^L^{r, x)^g (r, x)+ 5 (r, a;) + • • • , 



(40) 



(41) 



(0) (1) 

(p{r,x)—(p (r, x)+ (p {r,x) + 



We impose the boundary conditions at the position of the positive tension brane 



(n) 

g {r — 0, x) — 0, for all n > 0, 



(42) 



(43) 



(«) 

(p (r = 0, x) = 0, for all n > 0. (44) 



Other quantities are naturally expanded as 



(0) (1) 

Kf,^ =kf,y + kf,y +■■■ . (45) 
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2. 0* Order (Background geometry) 

The 0*'* order system can be easily integrated with respect to the extra dimension coordinate r to get the particular 
solution 



(0) (0) 



(46) 



(0) (0) 

^ (r,x) =if> {x). 



(47) 



This solution clearly satisfies the order junction conditions 



3. 1^* Order 



At 1** order the evolution equations are 



1 ^ (0) _962 , 1 / OJ,2 _ C (1) 1 ^ (0) „6l,2+2 ^/9/) 1 ^ ~9b^+l (1) 

[sb^ + i " ab^ + i " r I 3b^ + i ^' " 



-H 362+1 ^ ^ ^ 

' I" 3&2 + 1 I" 



(0) (0) 

V2& ^j^ + (262 + 1) ^iM,^!^ 



V26 -3f-3 



362 + r 



/(") (0) \ (0) 

(36^ + 1) cp\^ + ^1" q„<5^: 



+ 7? .9 =0 



(0) 



(48) 



1 362 - 3 
7362 + l' 



(0) 



-2^26 ip 



-1 (y 



(0) 



H 36^+1 X +e" 



-2V26 



9 (1) 

-66*^+2 1 ' 



1 2V2b 2V2b 



H 3i.2 + l 



/ 362 - 1 



C, 



Z 362 + 1 

(0) 



:?) -9bS+i (1) 



362 + 1 I" 362 



fi^/263 \ 2 / (0) 



(49) 



-2v2fc ^ 36^^,1 



(0) 



jl) 1 1 



(1) 



(1) 



+7 352 3\/26 if + (962 _ 5^ 



(50) 



The junction conditions at this order are 



■(1) 



r=0+ 



(51) 



■(1) 



r=0+ 



(52) 



In the preceding equations all the indices are raised with the zeroth order metric. Combining the trace of equation 
(I48|l with equation 149|1 we obtain 



1 ^ (0) (1) 



362 + 1 
6 



, (0) 

R\ 9 



(0) (0)' 



^24^ V26 _ 

iJ 3b2 + i _g 



^ (0) (1) 

<^^, . (53) 
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Imposing the junction conditions H51I52|I we get 



and the equation of motion for the 4D effective scalar field 



'(o)\ (0) (0) 
Rig) 



(0) 



Equation 153|) now reads 



Using the decomposition of Ki^ in 



equation 148|l can be easily integrated to find 



cl2 + V2&q, ^i"= 0. 



^ (1) 

-V26 (f^r ■ 



(1) (1) 



1 



(1) 



(54) 



(55) 



(56) 



(57) 



-2\/2b 



(1) 



(0) 



(0) 



(0) 



(0) 



(0) (0) 



traceless 



(58) 

where x'ii^) is an integration constant and the subscript [ ]traceiess means the traceless part of the quantity between 
square brackets. In terms of , the junction conditions (|51|l are 



r=0+ 

r=r 



0. 



(59) 



From the previous junction conditions H59(l we can obtain the 4D effective equations of motion 



R'^Ch = [C-- 



(0) 



1 / \"/ \"/ 



W 4 '' la 



(0)\1 /(O) , (0)' 



(0) (0) l2 (0) (0) 

(262 + 1) ^i^^h 



and 



(60) 



(61) 



Combining equation H56|) with the scalar field equation H5U|) and integrating it once with respect to the extra 
dimension, we get (after using the previous equations of motion to simplify the result) 



O AtJ?' -3I,>5 (1) 



r 



(0) 



(0) 



(0) \ (0) 



~p\l +a/26^|„^I" -^l"q,r + S(x), 



(62) 



where ,:i(a:) is just an integration constant. 
The junction conditions (|52|1 give 

S(x) = 

and the equation of motion for the second 4D effective scalar field 



(0) (0) (0) / 1 \ -1 (0) 



(63) 



(64) 
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C. The 4D Effective theory 

The 4D effective equations of motion are summarized as 



1 , / (0) (o)^ 



(0) 



(0) 



(0) 



(0) 



+ (262 + 1) ^i^^lA'.-j^; 



, (0) . (0) 



(65) 



(0) 



c\l + V2&q„ 0, 



(66) 



(0) (0) (0) 



-1 (0) 

^1" C, 



and they can be deduced from the foUowing action 



(0) , 

5 



(0) 



2'^ 



V2fc 



(0) (0) (0)' 

ii;( )- 



(67) 



(68) 



where la denotes covariant derivative with respect to 5^^. We should note that this effective action can be derived by 
substituting in ||SJ) the 5D solutions up to the first order and integrating it over the fifth dimension 
As a consistency check, we see that if we perform the conformal transformation 



(0) 



h^,^{x) = C^^(a;) ~g^,y (x), 
the previous action reduces to (|13|l if the effective scalar fields of the two theories are related through 



= C 



(69) 



(70) 



(71) 



The check consists in seeing that these relations are exactly the ones required so that the two observables (the proper 
distance between branes and the scalar field on the positive tension brane) agree in both approaches. 



D. The 5D exact solution 



It is straightforward to find a cosmological solution of this 4D effective theory. For example we can easily find the 
following particular solution 



(0) (0) 

giiu {x) = rif,^, C{x)=ht, (x) = 0, 



(72) 



where h is an integration constant. 

Now we are ready to address the question why the above solution can be lifted to an exact 5D solution. Let us start 
by calculating the next order correction '^^^g^v and ^^'^(p. Eq. (|62() and the boundary conditions 144(1 give '^^-'i^ = 0, if 
we take as 0*'' order solution Eqs. (|72|l . We can construct '^'^^K^n from Eqs. H56|l and (|58|l . For the 0*'* order solution 
(f7^ this gives ^^'K^u — 0. After imposing the boundary conditions ((43|l . we obtain that the next order correction 
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vanishes, ^^'^gfj,v{r,x) = 0. For solution (|72() it turns out that all the corrections vanish and the 0*^ order solution is 
an exact solution of the non-perturbed 5D Eqs. H34I39|I . 

For other solutions of the 4D effective theory, higher order corrections will not vanish and therefore they should be 
taken into account in the reconstruction of the 5D metric. Using the gradient expansion method, we can reconstruct 
the 5D solution pcrturbatively. We should emphasize that the choice of the 0*'* order metric is quite important in 
order to reconstruct 5D solutions efficiently. Our metric ansatz has the advantage that it is possible to recover the 
exact solution of Chen et al. (|72|l at 0*'* order. Indeed, if we had started with an ansatz like Eqs. H11I12|) we would 
need an infinite number of higher order terms to obtain the exact 5D solution. 

V. VALIDITY OF 4D EFFECTIVE THEORY 

In this section we will make comments on a recent work by Kodama and Uzawa |13|. Let us start by briefly 
describing their arguments. After deriving the 4D effective theory for warped compactification of the 5D Hofava- 
Witten model (they also extend their analysis to lOD IIB supergravity and obtain the same conclusions), the authors 
show that the 4D effective theory allows a wider class of solutions than the fundamental higher dimensional theory. 
Therefore we should be careful in using this effective theory approach, because we may find 4D solutions that do not 
satisfy the equations of motion once lifted back to 5D. 

The authors assume a metric ansatz of the form 

dSl = h{x,r)dr^ + {x,r)g^,^{x)dx^dx'' , (73) 

where the warp factor has the form r) = C{x) — r/L. This corresponds to taking K^^i, = 0, = and 6 = 1 (for 
the Hofava-Witten case) in our work. Then Eq. H34|l reduces to 

-il(r,x)-i (C{x)\: + \c{x)\l5'i^ + Hir,xr^ R^^gix)) = 0. (74) 

In order to satisfy this equation for all values of r, we should have 

R^^ig) = 0, q^, = 0. (75) 
They obtain the 4D effective action, by integrating the fifth dimension, as 

Seff (X I d^x^g {c{x) ~ i?(5), (76) 

which agrees with our effective action (I68II . As we have shown, this theory admits solutions with R.f^i,{g) ^ (see 
Eq. I|t)5|) '). which do not obey the constraint (|75|l obtained from the 5D equations of motion. 

However, it is clear from our analysis that their metric ansatz is too restrictive. If we consider a more general metric 
as our metric ansatz, we see that the 5D Einstein equations contain more terms given by K^^. With the inclusion 
of these new terms, the 5D equations do not necessarily imply H75|l . Of course, the non- vanishing if^^ changes the 
metric (|73|) and one could argue that the resultant 4D effective action would be also changed. However, it is shown 
that even if we include the first order corrections ^^^gfiv(x) to the metric, the resultant 4D effective action derived 
by integrating out the fifth dimension does not change [Tj. Therefore, for 4D solutions that do not satisfy H75|) . we 
should include the corrections to the metric ifT^ . We have provided this correction pcrturbatively. Using Eqs. 
and (|58|l . we can reconstruct the correction to the metric, *^^^g^i/(x), which is necessary to satisfy the 5D equation of 
motion. 

We should emphasize that the validity of the 4D effective theory is based on the conditions (|40|) . If the 4D effective 
theory admits a solution that violates the conditions (|40|) . then there is no guarantee that the 4D solution can be 
lifted up to the 5D solution consistently. We should check the validity of the 4D effective theory by calculating the 
higher order corrections to ensure that the higher order corrections can be neglected consistently. 

VI. CONCLUSION 

In this paper, we studied the moduli instability in a two brane model with a bulk scalar field recently found by 
Chen et al.. This model can be viewed as a generalization of the Hofava-Witten theory and the Randall-Sundrum 
model. The scalar field potentials in the bulk and on the branes are tuned in order to satisfy the BPS condition. 
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We used a low energy effective theory, which is derived by assuming that variations along the brane coordinates of 
the metric are small compared with variations along the dimension perpendicular to the brane. The effective theory 
is a bi-scalar tensor theory where one of the scalar fields arises from the bulk scalar field (dilaton) and the other arises 
from the degree of freedom of the distance between branes (radion). In the Einstein frame, the theory consists of two 
massless scalar fields, and the lack of potentials for these moduli fields was shown to be responsible for the instability. 

We found that the exact solution derived in 6J can be reproduced from the 0*'' order of the perturbative method, 
despite the fact that slow-motion approximations are used. We revisited the gradient expansion method which is used 
to derive the effective theory, in order to understand why the exact solution derived in ^ can be reproduced within 
the 4D effective theory. We proposed a new metric ansatz which is useful to see the relation between the solutions in 
the effective theory and the full solutions for 5D equation of motion. Using this metric ansatz, it is transparent why 
the moduli instability solution can be lifted to a full 5D solution. We have also shown that not all solutions in the 4D 
effective theory can be lifted to exact 5D solutions. For these solutions, the solutions in the effective theory receive 
higher order corrections in velocities of the branes and we need to find 5D solutions perturbatively. 

Finally, we comment on the recent arguments against the 4D effective theory. Ref. ^3] claims that the 4D effective 
theory allows a much wider class of solutions than the 5D theory. We argued that this conclusion comes from a too 
restricted metric ansatz used in Ref. 13]. Using a more general metric ansatz, we provided a way to reconstruct the 
full 5D solutions from the solutions in the 4D effective theory. 

Our method can be applied to other warped compactifications such as lOD type IIB supergravity models. In fact, 
there have been debates on the validity of the metric ansatz commonly used to derive the 4D effective theory. Our 
lOD generalization of the 0*'* order metric ansatz agrees with that proposed in Ref. |0| and the method presented 
in this paper will provide a consistent way to reduce the lOD theory to the 4D effective theory based on this metric 
ansatz. This will be reported in a future publication [l5|. 
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